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We consider the Ising model and the directed walk on two-dimensional layered lattices and show 
that the two problems are inherently related: The zero-field thermodynamical properties of the 
Ising model are contained in the spectrum of the transfer matrix of the directed walk. The critical 
properties of the two models are connected to the scaling behavior of the eigenvalue spectrum of 
the transfer matrix which is studied exactly through renormalization for different self-similar dis- 
tributions of the couplings. The models show very rich bulk and surface critical behaviors with 
nonuniversal critical exponents, coupling-dependent anisotropic scaling, first-order surface transi- 
tion, and stretched exponential critical correlations. It is shown that all the nonuniversal critical 
exponents obtained for the aperiodic Ising models satisfy scaling relations and can be expressed as 
functions of varying surface magnetic exponents. 



I. INTRODUCTION 

The study of layered Ising models (IM's) has been an 
active field of research during the last decades. One may 
mention the, pioneering workS|-pn two-dimensional (2D) 
periodicalljU and randomlyoQ layered lattices. Sim- 
ilarly the critical j-behavior of directed walks (DW's) 
in inhomogeneouil or random mediaQ has attracted 
widespread interest. 

Recently, following the discovery of quasicrystals,Ll on 
the one hand, and the progress in molecular beam epitaxy 
which allows the production of good quality multilayers, 
on the other handjj there has been a growing interest in 
the theoretical study of phase transitions in quasiperU 
odic systems and, more generally, aperiodic systems.B 
These are deterministic but nonperiodic structures which 
are called quasiperiodic when the spatial fluctuations are 
so weak that the Fourier spectrum is still discrete, but 
point symmetry is incompatible with a periodic struc- 
ture. Such systems may be considered as intermedi- 
ates between homogeneous and random ones and, conse- 
quently, are expected to display a rich variety of critical 
behaviors. 



A. Previously known results 

Most of the early works about phase transitions on 
aperiodic systems were done on quasiperiodic lattices 
and did not show any sign of modified critical behavior. 



Among these, one may mention an approximate renor- 
malization grcpp treatment of the classical IM on the 
Penrose latticald and MontCpGaiJo renormalization group 
studies of the same problemllil.ll3 Universal behavior was 
also obtained in Monte Carlo simulations of the peijco- 
lation problem on the Penrose lattice and ii|a-dualE3 as 
well as for the statistics of self-avoiding walks.Lj One may 
notice, as an exception, the analytical renormalization 
group study of interfacial roughness, still on the Pen- 
rose lattice, where the fluctuating interface "feels" a Fi- 
bonacci quasiperiodic potential. In this case, a marginal 
behavior was-, obtained for the decay of the transverse 
correlations. ta 

Probably the most studied system is the aperiodically 
layered 2D classical IM and j±s quantum counterpart in 
the extreme anisotropic limitjl3 the aperiodic Ising quan- 
tum chain in a transverse field. 

In the classical formulation, the energy of a configura- 
tion is given by 

<7k,icrk+l,l , 

(1.1) 

where the tr's are the spin-1/2 Ising variables, and Ki 
and K2 are the exchange interactions in the vertical and 
horizontal directions, respectively. Their values are the 
same in a vertical layer k and are modulated according 
to some aperiodic sequence in the horizontal direction. 

In the extreme anisotropic limit (Ki — > 00, K2 — *■ 0), 
the transfer matrix between successive rows in the verti- 
cal direction can be written as exp(— rTi), where t — 2K* 
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is the infinitesimal lattice spacing in the Euclidian time 
direction. TL is the Hamiltonian of a spin-1/2 quantum 
Ising chain in a transverse field: 
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L 

k=l 



L-1 



hkCTl-lY.-^'^^k^k+i, (1-2) 



fe=i 



where the c^'^'s are Pauli spin matrices. The transverse 
field hk (such that hkKl is the dual coupling Klik) given 
by eyii>[-2Kl{k)] = tanhii'i(fc)) plays the role of the 
temperature. The coupling is the ratio K2{k)/Kl. 
In general, due to universality, the classical and quan- 
tum systems have the same critical properties. Only in 
the presence of a marginal perturbation which leads to 
nonuniversal exponents are the corresponding quantum 
exponents obtained by taking an appropriate limit of the 
classical expressions. 

Usually the aperiodic modulation is assumed to involve 
only the horizontal couplings, i.e., the two-spin interac- 
tion Jfc and /ifc = 1. The fluctuation around the average 
J at a length scale L is measured by 



A(L) = ^(J, - J) 



(1.3) 



k=l 



When the aperiodic couplings are generated via substi- 
tutions using an inflation rule, this quantity behaves as 
L", where w is the wandering exponent of the aperiodic 
sequence wich is linked f.n.the two leading eigenvalues of 
a substitution matrix.EZllla 

For sequences with bounded fluctuations (w < 0) the 
aperiodic perturbation does not change the Ising criti- 
cal behavior. This was shown analytically by Tracy in 
the case of the Fibonacci sequence, with a; = — 1, for 
the 2D layered IM. The Onsager logarithmic singular- 
ity of the saecific heat then keeps a nonvanishing finite 
amplitude^ The same conclusion was reached for the 
Ising quantum chain with generalized Fibonacci modula- 
tions of the couplings. E] The low-energy spectrum of the 
quantum chain which, through the gap-exponent relation 
on finite critical chains, gives the values of the critical 
exponents was sliown to be unaffected by a quasiperi- 
odic modulationEd Universal behavior was also obtained, 
with the Thue-Morse sequence and its generalizations.c3 
In this case, the quantum chain is not quasiperiodic but 
the fluctations remain bounded. 

For an aperiojdic sequence with unbounded fluctuations 
(w > 0), TracynJ noticed that the Onsager singularity 
is suppressed like in the randomly layered McCoy- Wu 
model.H _ 

The situation was later clarified by LuckEj who 
proposed a generalizalion of the Harris criterion for 
quenched randomnessEj adapted to the case of aperiodic 
fluctuations of the couplings (see also Ref. 26 ) . By com- 
paring the mean shift of the local temperature in the 2D 
layered system (governed by the wandering exponent lo), 
at the scale of the correlation length of the unperturbed 



system, to the deviation from the critical temperature, 
one obtains a crossover exponent (j) = l + viu)—!). It 
controls the evolution of the amplitude of the aperiodic 
modulation when one approaches the critical point. For 
the 2D IM with v =1, the crossover exponent is equal 
to the wandering exponent so that, quite generally, the 
aperiodic modulation becomes a relevant perturbation 
and changes the Ising critical behavior when the flucttu. 
ations are unbounded, as conjectured by Benza et a/.cJ 
One must notice that the correspondence between rel- 
evant perturbations and unbounded fluctuations holds 
only when v = 1. The marginal behavior obtained for 
the Fibonacci sequence with bounded fluctpations in the 
case of the interface roughness problemlla follows from 
Luck's criterion where v is now the correlation length 
exponent in the transverse direction = 1/2. Finally 
let us mention that for a randomly layered system the 
relevance-irrelevance criterion applies with a; = l/2. 

In the same workj2J Luck checked the validity of 
his critierion for random and aperiodic quantum Ising 
chains. To treat the aperiodic problem he considers 
periodic approximants, i.e., a periodic quantum chain 
with a large unit cell of length L in which the couplings 
Jk are distributed according to the aperiodic sequence. 
He deducejS- the low-energy behavior of the fermionic 
excitationsEZI A from a perturbation expansion in A. For 
the unperturbed problem at criticality, the massless exci- 
tations have a linear low-energy dispersion relation A~vq 
where v is the velocity and q the wave vector. On the pe- 
riodic approximant a L-dependent velocity vl is obtained 
and the properties of the aperiodic system are governed 
by the limiting behavior of when L goes to infinity. 
The behavior of the singular part of the ground-state en- 
ergy, corresponding to the free energy in the 2D classical 
system, is linked to the low-energy excitation spectrum 
and its temperature dependence can be obtained through 
a scaling argument. 

For sequences with bounded fluctuations (a; < 0), vl 
is bounded and nonvanishing in the limit L — s- oo so that 
the Onsager logarithmic critical singularity is preserved. 
For unbounded fluctuations (a;>0), the typical velocity 
vanishes exponentially, leading to an essential singularity 
for the singular part of the ground-state energy as in the 
case of random chains. Finally, when the fluctuations 
grow on a logarithmic scale {lu — Q), the typical velocity 
vanishes as a nonuniversal power of L. The perturbation 
is marginal and the specific heat exponent is negative (the 
logarithmic singularity is suppressed) and varies contin- 
uously with the amplitude of the aperiodic modulation. 
This marginal behavior was checked numerically. 



B. Renormalization-group method and main results 

The results obtained so far for different aperiodic mod- 
ulations in diffofertt-, models are all in accordance with 
Luck's criterionJ23n2a. Most of the activity in our groups 
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was concerned with the study of the surface and bulk 
critical properties of 2D aperiodically layered IM's, ei- 
ther using the ID quantum formulation or working on 
a triangular lattice, making use of the star-triangle rela- 
tion. 

Although some relevant perturbations were treated in 
Refs. p8| and |29|, we mainly considered marginal aperiodic 
perturbations. The continuously varying surface mag- 
netic exponent Xm^ = f^s/i^ was obtained analytically for 
different aperiodic sequences whereas the scaling dimen- 
sion of the surface energy waSr-Gpnjectured on the basis 
of finite-size scaling studies.E3~Ea 

The marginal aperiodic-|Hipdels were found to dis- 
play anisotropic scalingj23'EllE3 The correlation length di- 
verges with different exponents along and through the 
layers with a ratio z — giving a continuously vary- 

ing anisotropy exponent. Such a jbehavior was in fact 
implicitly contained in Luck's workcJ where a power-law 
dependence on L was found in the marginal case. Accu- 
rate numerical calculations of the anisotropy exponent z 
led us to propose a simple scaling relation between z and 
the surface magnetic exponents on both sides of the sys- 
tem, Xm, and Xm, ■ The anisotropic scaling of bulk and 
surface properties was extensively studied in Ref. 

In this paper we present the results of an exact 
renormalization-group (RG) study of aperiodic and hier- 
archical Ising and DW models. The introduction of RG 
techniques into the field of phase transitions and critical 
phenomena has largely contributed to our understand- 
ing of the properties of the critical state. For instance, 
the RG method has given a natural explanation for the 
scaling hypothesis and universality. At the same time, 
it has provided, powerful procedures to calculate criti- 
cal exponents,E3 albeit generally using some approxima- 
tion, e.g., approximate RG transformations, expansions 
in a small parameter (e, 1/N,...), or numerical meth- 
ods. There are few nontrivial problems in statistical me- 
chanics for which the RG transformation can be worked 
out exactlji One may mention the IM on the triangu- 
lar latticeEZl or different physical processes on self-similar 
fractal objects.L3 

Here we develop exact RG solutions for a class of 2D 
layered Ising and DW models. The novel feature of 
our approach is that we study both problems within the 
framework of the same RG transformation. It is based 
on a hitherto unnoticed connection between the eigen- 
value problem for fermionic excitations which enters the 
solution of the IM (Ref. |27|) and the transfer matrix of a 
DW in two dimensions. Both problems are considered on 
layered lattices, such that the walk is directed along the 
translationally invariant direction. The solution of the 
DW, which means the diagonalization of its transfer ma- 
trix (TM), provides in principle all the necessary informa- 
tion to obtain the zero-field thermodynamical properties 
and correlation functions of the IM. 

The critical properties of the two models are connected 
to the scaling behavior of the eigenstates of the TM at 
different edges of the spectrum. An exact RG study of 



the eigenvalue problem of the TM is performed for dif- 
ferent self-similar distributions of the couplings and the 
critical properties of the IM and the DW are governed by 
two different fixed points of the same RG transformation. 

Our method is well adapted to the case of self-similar 
perturbations. It is quite different from the approximate 
renormalization group technique recently introduced by 
Fisher to treat randomly layered systems In this ap- 
proach, which leads to exact results in the critical do- 
main, instead of using the transformation to fermions, 
Fischer works on the Hamiltonian itself, reducing the en- 
ergy scale of the problem by a systematic elimination of 
the stronger couplings 

The bulk and surface critical properties are examined 
for several aperiodic and hierarchical sequences. In all 
of the models we studied the aperiodicity is marginal at 
the Ising fixed point and induces continuously varying 
critical exponents. The bulk anisotropy exponent z and 
the correlation length exponents v, on the one hand, and 
the surface energy exponent Xe^ and the surface mag- 
netization exponents Xm^ and Xm^^ on the other hand, 
are obtained analytically. We also prove the previously 
conjectured relation between the anisotropy and surface 
magnetic exponents, z = Xms+Xm, > 1, which holds for 
sequences which modify the critical coupling. A simple 
scaling picture emerges in which the surface magnetic 
exponents play a fundamental role: All the nonuniversal 
exponents (except the bulk magnetic one, not considered 
in this work) can be expressed as functions of these two 
surface magnetic exponents. The surface energy expo- 
nent is given by Xe^ = z + 2xm^ on one side of the chain 
and Xe^ = z+2xm^ on the other, whereas the specific heat 
exponent is given by a = 1 — z < 0. 

With aperiodic sequences with a vanishing density of 
modified couplings, which changes the critical behavior 
only locally near the surface, there is no anisotropy in the 
bulk of the system (z = 1). For sufficiently strong mod- 
ified couplings, the surface remains ordered at the bulk 
critical point and then the first excitation alone scales in 
an anomalous way, with a continuously varying power of 
the size of the system. 

Finally one may mention that these marginal aperi- 
odic Ising sytems are closer to periodic than to randomly 
layered ones. The varying exponents evolve continuously 
from their unperturbed values when the aperiodic mod- 
ulation grows. In particular we checked numerically that 
the gap is nonvanishing in the disorderediihase, i.e., that 
there is no trace of a Griffiths phascxja as expected 
for systems with bounded fluctuations .c3 Even with rele- 
vant aperiodic perturbations, which by some aspects are 
closer to random ones, displaying essential sirigularities 
in the singular part of the-ground-state energycJ and in 
the surface magnetization, ESLthe Griffiths phase is absent 
according to a recent study.Cj 

A short account of our results, concerning thCj-hulk 
critical behavior, has been given in a recent Letter 

The structure of the paper is the following. The rela- 
tion between the IM and the DW is presented in Sec. ||, 
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afterwards the basic pro per ties of aperiodic sequences 
are reca pitula ted in Sec. [II. In the foUowing sections, 
Sees. [V VIII, the RG transformation is worked out for 



different aperiodic and hierarchical models which were 
chosen in order to illustrate the different renormalization 
procedures one may use, in particular in the treatment 
of the surface properties. Both the bulk and surface crit- 
ical behaviors are studied. Some relations between the 



critical exponents of the IM are derived in Sec. IX and 
the results are discussed in Sec. ^ Details about the 
derivation of the RG equations for the specific models 
are collected in the Appendixes. 



canonical transformationc3 which gives 



9=1 



(2.1) 



where rj^ and rjq are fermion creation and anihilation op- 
erators, respectively. The fermion excitations are non- 
negative and satisfy the set of equations 



-Jfc-l*,(fc-l)-/lfe*g(fc) 



(2.2) 



II. RELATION BETWEEN THE ISING 
QUANTUM CHAIN AND THE DIRECTED 
WALK MODEL 

Using a Jordan- Wigner transformation,^! the Ising 
Hamiltonian (l.S) can be rewritten as a quadratic form 
in fermion operators. It is then diagonalized through a 



with the boundary conditions Jq — Jl — 0. The vec- 
tors $g's and ^q's, which are related to the coefficients 
of the canonical transformation, are normalized. They 
enter into the expressions ot-iCorrelation functions and 
thermodynamical quantities .E3 

Usually one proceeds by eliminating either 'i'q or $q 



in Eqs. ( |2.2| ) and the excitations are deduced from the 
solution of one of the following eigenvalue problems: 



Jk-ihk'Hqik- 



- 1) + i4_, + hi)^,{k) + Jkhk^ik - 

1) + (Jk + hl)-^q{k) + ,hhk+i-^q{k - 



^l%{k) 
A^*,(fc) 



(2.3) 



with the same boundary conditions as above. This last 
step can be avoided by introducing a 2L-dimensional vec- 
tor Vq with components 

t/,(2fc-l) = -<f,(fc), V,{2k)^-^,{k) , (2.4) 



and noticing that the relations in Eqs. ( p. 2D then corre- 
spond to the eigenvalue problem for the matrix: 



/O/iiO 0---\ 
/ii Ji ••■ 



Ji /i2 
/i2 J2 



v 



(2.5) 



Taking the square of T, odd and even components of Vg 
are decoupled, and one recovers the two eigenvalue equa- 
tions in Eqs. ( ^.3|) . The matrix T can be interpreted as 
the TM of a DW problem on two interpenetrating, di- 
agonally layered square lattices as shown in Fig. |l|. The 
walker makes steps with weights and Jk between first- 
neighbor sites on one of the two square lattices and the 
walk is directed in th e di agonal direction. 

According to Eqs. (2^), changing into — $g in V^, 
the eigenvector corresponding to — Ag is obtained. Thus 
all the information about the DW and the IM is con- 
tained in that part of the spectrum with Ag > 0. Later 
on we shall restrict ourselves to this sector. 




FIG. 1. Portion of two interpenetrating diagonally layered 
square lattices involved in the transfer matrix of the directed 
walk. The fugacities are alternatively and Jk- 

Let us now consider the correlation lengths in the di- 
rection parallel to the layers in both problems. For the 
DW it can be expressed as a function of the two leading 
eigenvalues of the TM with 



In 



A, 



A 



A, 



Ar,-A 



L-l 



(2.6) 



Thus ^1*^ is proportional to the inverse gap at the top 
of the spectrum of the TM. For the IM in the disordered 
phase, the correlation length is the inverse of t he lo west 
excitation energy of the Hamiltonian Ti. in Eq. (2.1): 



IM 



Ar 



The fermionic excitation Ai is al| 
in the positive sector of the TM.I 



(2.7) 

the lowest eigenvalue 



4 



When any one of the two critical points is approached, 
the correlation length of the problem diverges and the 
corresponding part in the TM spectrum displays a scal- 
ing behavior. Let us consider a finite system with trans- 
verse size L^l and denote by A A; either A^ — Ai_i for 
the DW or A^ for the IM with i-^L. Under a change of 
the length scale by a factor b>l such that L' = L/b, the 
gaps are expected to transform as 



(AA,)' = b^^ AA, 



(2.8) 



where the scaling dimension is generally different for dif- 
ferent parts of the spectrum. This leads to the finite- 
size scaling be havio r AAi(L) ^ L^^^ and, according to 
Eqs. ( |2.6|) and (2/7), the longitudinal correlation lengths 
behave as ^|| ^ L^^. Since the anisotropy expo- 

^G0i: 



nent z, defined through ^|| 



is given by 



I/A 



(2.9) 



In the case of the DW one is interested in the transverse 
fiuctuations of the walk which are characterized by a wan- 
dering exponent w through ~ f ™ (Refs. ^ and so 
that 



(2.10) 



III. SUBSTITUTION MATRIX AND 
RELEVANCE-IRRELEVANCE CRITERION 

In the following, except for the hierarchical sequence, 
we consider sequences generated via substitutions on a 
finite alphabet such that, in the case of two letters A and 
B, A^S{A), B^S{B). The propertiejs-|o£,the sequence 
are governed by its substitution matrixtlrEJ 



M 



where the matrix element nf gives the number of i in 
S{j). The matrix elements in M" give the same num- 
bers in the sequence obtained after n iterations. When 
the substitution starts with j, the corresponding numbers 
are contained in column j. 

If Ujy denotes the right eigenvector of M with eigen- 
value f2,y, the asymptotic density of i is given by 



(3.1) 



(3.2) 



where Ui is the eigenvector corresponding to the leading 
eigenvalue fii. The length of the sequence after n itera- 
tions is related to the leading eigenvalue through L„ ^ il" 
so that fii > 1. 

In the following, each letter in the sequence is replaced 
by one digit or more (for examples, A — 0, B = l). Thus 



one obtains a sequence of digits fk (fc = l,2, . . . , L). The 



aperiodic Hamiltonian is defined as in Eq. (1.2) with a 
constant transverse field hk — h and a modulation of the 
couplings following the aperiodic sequence. 



(3.3) 



where J is the unperturbed interaction and R the mod- 
ulation ratio. 

When fk — 0, 1, the cumulated deviation A(L) from 
the averaged coupling J defined in Eq. ( p_.3| ) scales with 
L as 



A(L) = J(i? ~ l){nL - Lpoo) ~ 5\n2\ 
In this expression. 



L 

k=l 



riL 

= lim — 

L — >oo L 



(3.4) 



(3.5) 



give the number of digits equal to 1 in a sequence with 
length L and their asymptotic density, which can be de- 
duced from Eq. ( |3.2| ), respectively; ^2 is the next-to- 
leading eigenvalue of the substitution matrix, 5 measures 
the amplitude of the aperiodic modulation, and lo is the 
wandering exponent of the sequence, given by 



InlOal 
Inf^i 



(3.6) 



Thus the mean shift of the coupling strength SJ{L) at a 
length scale L, proportional to L'^~^, is governed by the 
wandering exponent. 

The relevance of the perturbation follows when one 
compares the deviation t from the critical point to the 
averaged temperature shift 6t ^ 6J{^) induced by the 
aperiodicity a^-a length scale given by the correlation 
length ^~t-'':El 



5t 

t 



= 1 + J/(CJ - 1) 



(3.7) 



When 0>O, the ratio is divergent, which indicates a rele- 
vant perturbation. When (/)<0, the ratio vanishes at the 
critical point and the perturbation is irrelevant. Finally, 
when = 0, the perturbation is marginal and may lead 
to a nonuniversal behavior. The same conclusions can 
be reached by calculating the scaling dimcnsicp, of the 
modulation amplitude 5, which is equal to .E3 

Fo r a strongly anisotropic unperturbed system u in 
Eq. (^^) has to be replaced by the exponent of the 
correlation length in the-direction perpendicular to the 
layers in the 2D systems 

The critical transverse field of the inhomogeneous IM 
is generally given byca 



lim 



k=l 



(3., 
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Introducing the reduced coupling A = J//i, its criti- 
cal v alue o n th e aper iodi c quantum chain follows from 
Eqs. Q, O), and as 



A. 



(3.9) 



IV. PERIOD-DOUBLING SEQUENCE 

A. Definition and general properties 

The period-doubling sequenceEl follows from the sub- 
stitutions A S{A) = AB, B S{B) = AA. Here we 
make the identification A = and i? = 1, i.e., J a — J 
and Jb = JR according to Eq. (3.3). Thus starting on A, 



after n iterations, one obtains the following sequences of 
digits /fc: 

n = , 

n = 1 , 01 

n = 2, 0100 

n = 3, 01000101 

n = 4, 0100010101000100. (4.1) 

The eigenvalues of the substitution matrix are Sli = 2 
and = —1 so that the wandering exponent lo, given by 



Eq. ( |3.6D , vanishes. The asymptotic density p^—p^^o' 
1/3 follows from Eq. (3.2) and leads to the Ising critical 
couphng Ac = i?"^/^, according to Eq. (3.9). 

One easily verifies on the last line of Eqs. ( [4.l| ) that 
the fk's satisfy the relations 

/2fe = 1 - /fe , /2fe+i = . (4.2) 



B. Built critical behavior 

We now proceed to the exact renormalization o f th e 
eigenvalue equations, associated with the matrix (2.5), 
which follow from Eqs. ( |2.2| ) and (2^). We first treat 
the bulk problem on a semi-infinite system. To recover 
the period-doubling sequence of interactions after one 
renormalization step, we eliminate triplets of interactions 
(J, RJ, J) indicated by crosses in Fig. |[ 



RJ 

h J 
x|x x|x x|x 



XX X 



X XX 



FIG. 2. Matrix elements Tk,k+i as a function of k for the 
period-doubling sequence. Components of the eigenvector to 
be decimated out in the RG transformation are denoted by 
crosses. The heights of solid vertical bars indicate the strength 
of the couplings; the grey bars stand for the field. 



Using reduced coupHngs A^ = Jk/h and a reduced 
eigenvalue A = A/h, the RG equations, as derived in Ap- 
pendix are given by 



A' = A 



i?A3 



A' = 



i?A2 



(4.3) 



where c and d are defi ned in Eq s. (A2 ). 

According to Eqs. (A2) and (A4), the components of 
the eigenvectors transform as 



V'{2k) = V{8k) , y'(2fc + 1) = y(8fc- 
The bulk IM fixed point corresponds to 

A* ^ , A* = , 



which, using Eqs. (A2), leads to 

c* = -l, d* = 1 + R^/^ + R-^^^ 



(4.4) 
(4.5) 
(4.6) 



Thus the eigenvalues of the linearized transformation are 
given by 



OA' 
9A 



(i? 



1/3 



l,yt 



dX 



and, with 6 = 4, one obtains 

_ ln(j?i/3 + 
^" h^2 ' 



yt 



= 4, 
(4.7) 

(4.8) 



thus confirming the conjecture of Ref. ^ 

The top of the spectrum, which governs the behavior of 
the DW, scales to a fixed point with A^oo and A^oo. 
Thus it is convenient to write the RG equations in terms 
of the new variables k=:1/A, and a — A/\ leading to 



,R 



a = a [ 1 



>B 



A' A^ 

A^[C?- R^)[C? - 1)2 + 2At2a2 (1 _ „2) ^ ^4^2 ^ 

B = ia^- R^){a^ - 1) + {1 ~ 2a^) + . (4.9) 
From Eqs. (|4.9D the DW fixed point is given by 







i?2 



. 



(4.10) 



The scaling behavior at the DW fixed point is differ- 
ent for the homogeneous model, with i?=l, and for the 
aperiodic one, with R^l. First we start with the homo- 
geneous model, where the separatrix in the (A, k) plane 
is linear at the fixed point: A(k) = a*K. According to 
Eqs. (|4.9| ), a point with coordinates A = aK, k<^ 1, which 
lies close to the separatrix, is repelled by the fixed point 
as 



2a 



A' 



1 



1 



a2(a2 -2) 



, (4.11) 
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and thus a' = A'/k' is given by 

a' = a* - W + : 



(4.12) 



The fixed-point value a* — 2 determines the equation of 
the separatrix while the leading eigenvalue of the trans- 
formation, ei = da' /da\* = 16, is connected to the gap 
exponent through j/a = In ei/ ln6 = 2. Consequently, the 
wandering exponent is given by 



1 



w = — = 



(i? = 1) , 



(4.13) 



in agreement with known results. EI 

For the aperiodic model the separatrix has a quadratic 
dependence A(k) = when k<C 1, with /?* = {^/2R— 
2B?) / [1 — R^), in contrast to the linear behavior for the 
homogeneous model. The scaling behavior of a point 
with coordinates A = /?k^, k <C 1 , close to the separatrix, 
can be deduced from Eqs. (4.9) ^ and A' A, 

thus A'/k'~A^/^/k. Consequently, at a fixed k = k'<c1 
we obtain A'~ A^^^, which represents a strong repulsion. 
This type of scaling behavior is compatible with an es- 
sential singularity in the gaps at the top of the spectrum: 



AAi ~ exp(-CL'") 



(4.14) 



with a — 1/2 since the rescaling factor is & = 4. From 



Eq. (4.14) the parallel correlation length of the DW is 
given by (AAi)"! - exp(CLi/2)^ tj^^s ^he trans- 

verse fluctuations of the walk grow anomalously, on a 
logarithmic scale: 



[X(t)-X(0)]2)i/2_ln2(t) 



(4.15) 



Here X(t) denotes the position of the walker at time t. 
We note that the same asymptotic behavior is found in 
the Sinai model of a |e*ie-dimensional random walk in a 
random environment e3 



C. Surface critical behavior 

We now turn to the renormalization of the surface 
block, looking for the scaling behavior of the surface 
temperature ts- In order to do so, we apply a modi- 
fied transverse field hi = hts on the first site. As shown 
in Appendix ^ the RG transformation now generates an 
auxiliary variable ^, in terms of which the recursion re- 
lations are given by 



c-d 



:_d_ 



(4.16) 



where c and d are the parameters defined in Eqs. (A2). 
The auxiliary variable 9, which does not enter into the 
renormalization of ts, may be discarded. 

For the bulk Ising fixed point values of c and d given in 
Eqs. ( |4.6| ) , the transformation of the surface temperature 
gives two surface fixed points with: 



dtf 
dtl 
dtf 



= (i?i/3 
= (i?i/3 



i?-l/3)-2 < 1 , 



.*2 



R 



)^>1, 



,*2 



(4.17a) 
(4.17b) 



Thus the attractive fixed p oint in the critical surface, 
corresponding to Eq. ( 4.17a ), governs the surface critical 
behavior and leads to the scaling dimension of the surface 
temperature, 



ln{R^/^ + i?-i/3) 



In 2 



(4.18) 



in agreement with the conjecture of Ref. ^ The same 
quan tity a t the repulsive fixed point, corresponding to 
Eq. (4.171:), is given by 



2 In 2 



(4.19) 



We now consider the critical behavior of the surface 
magnetization of the aperiodic IM. For a semi-infinite 
layered system, the surface magnetization is simply given 
by the first component of the normalized eigenvector #i 
corresponding to the lowest femiionic excitation, which 
vanishes in the ordered phase:La 



= $i(l), ^<i>2(fc) = l 



(4.20) 



k=l 



According to Eq. (2.4), its scaling dimension can be 
deduced from the renormalization of the odd components 
of V at the Ising fixed point which corresponds to A* = 0. 



Like V{8k+7) in Eqs. ([A2[), the odd components inside 
each block can be expr esse d as functions of V{8k+1) and 
V{8k + 8) using Eqs. (|Al|b)-(|Allg). At the fixed point, 
taking Eqs. ( [4.4[ ) into account, one obtains 



$*(4fc+ 1) = $'*(A: + 1) , 
$*(4fc + 2) = $'*(fc + l) , 
$*(4fc + 3) = -R-^/^ $'*(/c + 1) 
$*(4fc + 4) = $'*(A: + 1) . 



(4.21) 



Thus the normahzation of leads to 

oo 4 cx 



j2J2'^*H^k+i) = (i?'/-'+i?-'/')'5]$'*'(fc+i) = 1 . 

(4.22) 



fe=0 1=1 



Near the critical point, the surface magnetization trans- 
forms as 



$'*(!) 



(4.23) 



7 



so that, using Eq. ( |4.22| ) 
with 5 = 4, 



TO, 



{R^/^ + R-^/^)ms and, 
ln(i?iA^ + i?-i/3) 



21n2 



(4.24) 



in agreement |5dth an analytical result for the surface 
magnetization.Bj The thermal and magnetic surface scal- 
ing dimensions are related through 



1 



:yts 



(4.25) 



a relation conjectured in Ref. |3^. Furthermore, compar- 
ing Eqs. (4.24) and (4.1£), one may verify that 



(4.26) 



These relations, which are generally valid for the IM, will 
be discussed in Sec. EX. 



V. HIERARCHICAL SEQUENCE 
A. Definition and general properties 



One may notice that the underlined terms f2k give 
back the original sequence with /j; replaced by /^-l-l. 
The same property remains true for fmk with any value 
of TO so that: 



-1=0. 
(5.4) 



fmk fk ^" 1 7 fmk-]-l fnik+2 ' ' ' fmk^ 



Th e Is ing critical coupling is still given by Eqs. (|3.5| ) 
and ( ^.9[) where ul can be evaluated recursively. Using 
Eqs. (|5.4D for a sequence with L = mP, we obtain 



k=l 



,fe = rimP-i + TO' 



p-1 



rnF~^ -t- TO^"^ 



- 1 



which leads to 



TO — 1 



(5.5) 



(5.6) 



In the generalized hierarchical sequence associated 
with an integer m^J- the positions k of the digits fk 
satisfy the relationcBEll 

k = m-^'' (ml + p) , 1 = 0,1,... , p = 1,2, . . . ,m — 1 . 

(5.1) 

With TO 1^ 2, the Huberman-Kerszberg sequence is 
recovered.L2l 

We recently noticed that these hierarchical sequences 
can be also generated via substitution, using an alpha- 
bet with an infinite number of letters. Let us put the 
letters in correspondance with the natural numbers; the 
/fc's then follow from n-^S{n) withEil 



(5.2) 



m — 1 

S{n) = oFT^ (n + 1) 



Starting with n = 0, repeated applications of Eq. ( p.2\) for 
171 = 2 leads to the following sequence at the fourth step: 



0102010301020104 . 



(5.3) 



According to Eq. (3.3), it corresponds to the interactions 

J^ = J , Jo = J R, J?, = J , Ja = J R? , . . . . 



B. Bullc critical behavior 

In the exact renormalization group transformation, we 
decimate out those sites of the lattice, which have con- 
nection by a J coupling. In such a way, blocks of 2(to— 1) 
sites are eliminated as indicated by crosses in Fig. 0. Us- 
ing the reduced variables \ = J /h and A = A/h, a lengthy 
calculation detailed in Appendix ^ leads to the transfor- 
mation 



A' = A 



D 



2m-2 



D 



2m-3 



X^XR^^^, (5.7) 

A rrj, — 1 ' ^ ' 



where the D^s are determinants defined in Eq. (B3). 



RJ ^ 
h J 1 




|x|x| 1 |X|x 


x|x| 1 |x|x| 



FIG. 3. As in Fig. ^ for the hierarchical sequence with 
m = 2. 

With to = 2 and m = 3 one obtains 



A' = A (A2 -X'-l) 



y = R (A^ - A^) (to = 2) 



A' 



A^-A' 



2A2 + A^ + 1 



A' 



R 



2\2 



(A^ - A^) 



A^ 



(m = 3) 



(5.8a) 
(5.8b) 



The RG transformation in Eqs. (5.7) has an Ising fixed point with A* = 0. To study the behavior of the system 
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close to this fixed point, we expand the determinants -D2m-2 and -D2m-3 to linear order in A: 



D 



2m -2 



= (-A2)"-1 + 0(A2) 



D 



2m-3 



1 _ \2m-2 

= i-ir-^A——^ + o{p) 



Putting these expressions into Eqs. (5.7), we obtain the location of the Ising fixed point: 

A* = , A* = . 
The eigenvalues of the linearized transformation are given by 



9A' 

aA 



IA*| 



iA*r 



|A*|-|A*|- 



l,yt 



9y 

dX 



m . 



and with b = m one obtains 

_ ln(|i?™/(™-i) _ln(|i?i/(™-i)_i?-i/(™-i)|) 

Inm 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



in agreement with the conjecture of Ref. ^ 

Another fixed point of the transformation with A* > 
governs the critical behavior of the DW as shown in 
Fig. ^. The position of this fixed point and the values 
of the corresponding critical exponents can be calculated 
in a closed form only for m = 2 and m = 3. 



A 


DW 


IM 





FIG. 4. Schematic RG phase diagram for the hierarchical 
model. There are two nontrivial fixed points on the sepa- 
ratrix, at A* = for the IM and at a nonvanishing value 
of A* for the DW. The two fixed points generate fixed lines 
(not shown here), parametrized by the modulation ratio R, 
in an extended parameter space. These fixed lines govern 
the marginal behavior of the IM and the DW, respectively. 
For the period-doubling, three-folding, and paper-folding se- 
quences the scaling is anomalous at the DW fixed point which 
is shifted to infinity. 

T he D W fixed point for m = 2 is deduced from 
Eq. ( 5.8a ) as 



A* 



1 - R 



A* = 



R 



1 - R 



(5.13) 



The leading eigenvalue of the linearized transformation 
is given by 



R-^ +R+-] -2 



1/2 



and thus the wandering exponent of the walk is 
1 In 2 



UA In ei 



(5.14) 



(5.15) 



since the rescaling factor is 6 = 2. 

The DW fixed point for to = 3, which follows from 
Eq. ( 5.8b| ), is located at: 



A* 



A* = 



V2R 

l-R 



(5.16) 



l-R ' 
The leading eigenvalue reads 

ei = 2(i?-i +R+1)+ [4(i?-i +R+lf-3] , 

(5.17) 

and the wandering exponent of the walk is given by: 
1 In 3 



(5.18) 



VA Inei 
since the rescaling factor is now b — 3. 

C. Surface critical behavior 



At the surface of the system we define a modified sur- 
face field hi = htg with scaling dimension yt^ and in- 
troduce an auxiliary variable 6 to take into acount the 
asymmetry of the renormalized couplings. In terms of 
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the reduced variables, Eqs. (Bl) have to be supplemented 
for n = with the first two equations 



ts 



-AVil)+etsVi2) ^ , 

-F(l) - Ay(2) + AF(3) = . (5.19) 
Then, besides the recursion relations of the bulk variables 



in Eqs.( 5.7), we have two more relations for the surface 
fields: 

-,1/2 



ts 



D2m-2 




.^2m-2 - 


f t'^D2m-4 






D2m-2 4 


- tlD2m-4_ 



-,1/2 



(5.20) 



where Z?2m-4 is defined through D2m-3 — AZ?2m-4- 
Here, as before, the auxiliary variable 9 does not enter 
into the renormalization of tg and may be discarded. 

As one can see from Eqs. ( 5.20| ), there are two surface 
fixed points at t* =0 and at = 1, from which the latter 
is stable, both on the IM and the DW critical surfaces. 
Evaluating the linearized RG transformation around the 
stable fixed point, one obtains, for the scaling dimension 
of the surface temperature. 

In (1 + i?2/(™-l) + i?4/(m-l) + . . . + ^2) 



yt 



Inm 



(IM) 



(5.21) 

for the IM as expected from numerical resultsiil and 
InR 



Inm 



(DW) , 



(5.22) 



for the DW. 

The scaling dimension of the surface temperature for 
the DW is related to the anomalous diffusion expo- 
nent (i_|m the hierarchical lattice. According to exact 
results,EJ'L3 the mean-square displacement of a diffusive 
particle in a hierarchical environment is asymptotically 
given by {X'^{t)) ^t^^'^^ , where 



_ J 1 — In R/ him , 
2 , 



R < l/m 
R> l/m 



(5.23) 

for anomalous diffusion, i.e.. 



Thus one has = 1 + yt^ 
with R< l/m. 

One can also deduce the scaling dimension of the sur- 
face magnetization of the IM from the rescaling of the 
surface component of the eigenvector $i(l). Following 
the same way as for the period-doubling seque nce i n 
Sec. IV C, with the fixed-point parameters in Eqs. (5.10), 
one obtains 

ln(l + i?2/('™-l) + Ri/irn-l) + 
Xm, = ^-j (IM) , 

2 Inm 

(5.24) 

in agreement with Ref . ^ for m — 2 and Ref. ^ for any 
value of m. One can easily check that the scaling relations 



in Eqs. (4.25) and (4.26) are satisfied for the hierarchical 
IM too. 



VI. THREE-FOLDING SEQUENCE 



A. Definition and general properties 



The three-folding sequenceE^ is generated through the 
substitutions A S{A) = ABA and B S{B) = ABB. 
Starting on A with A = and B = l, at the third iteration, 
one obtains 



010011010010011011010011010. 



(6.1) 



The substitution matrix has eigenvalues fii = 3 and 
fl2 = 1, leading to the wandering exponent lu — 0. The 
asymptotic density is poc ~ 1/2 and gives Ac = i?^^/^ 
for the Ising critical coupling. The same sequence is re- 
covered when one keeps every third term, underlined in 
Eq. (^^). The digits in between are always and 1 so 
that the following relations are obtained: 



Hk — fk 



f: 



3fc+l 



= 



/: 



3k+2 



= 1 



(6.2) 



B. Bulk critical behavior 



To proceed to the bulk renormalization one considers 
blocks of six eigenvalue equations from which four, in- 
dicated by crosses in Fig. |^, are eliminated so that the 
rescaling factor is now 6 = 3. It is convenient to use the 
reduced eigenvalue A — A/ J and the temperaturelike pa- 
rameter II — h/ J as well as an auxiliary variable k which 
is needed to take into account the form of the couplings 
after renormalization. The decimation described in Ap- 
pendix ^ leads to the renormalized variables 



A' = A 
, Rf,^ 



(-D 



1-^ 

e/ \ e 



1/2 



12 2 

K = K ■ 



(6.3) 



where c, d, and e are the parameters defined in Eqs. ( |C2| ). 
The components of the eigenvector V transform as 



V'{2k) = V{6k) 



RJ 



y'(2fc+ 1) = y(6fc+ 1) . (6.4) 



h J 








X 1 X X 




X 1 X X 1 X 


X 1 X ' X 1 X 1 



FIG. 5. As in Fig. |2| for the three-folding sequence. 
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At the IM fixed point 

A* = , Ai* = , 
and, according to Eqs. (|C^), 

c* = -R^ (R+l), d* = -R{R+l), 



(6.5) 



e* = R^ 



(6.6) 



The auxihary parameter k, which does not enter into 
the renormalization of the physical variables, need not 
be further considered. The eigenvalues of the linearized 
RG transformation follow from Eqs. (|6.3|) with 



dA' 
9A 



[(2 + i?)(2 + i?-i)]i/2 , 



3 . (6.7) 



Thus, with 6 = 3, one obtains the bulk scaling dimensions 



\n[{2 + R){2 + R-^)] 
2 In 3 



yt = v 



= 1 



(6.8) 



as expected, according to the numerical study of Ref. |3^. 

The top of the spectrum again scales to infinity (A 
oo, J-^oo) such that the DW fixed point is located at 
A* = 1, /i* = 0. The equation of the separatrix is given 
by A = A^ — 1 = 7*/i when /i ^ 0, i.e., close to the fixed 
point. The scaling behavior of a point with coordinates 
A = 7/i, /i^l is of the form A'^A and like for 

the period-doubling sequence. Thus the highest gap in 
the TM spectrum displa ys an essential singularity of the 
form given above in Eq. ( 4.14 ) but (t = In 2/ In 3 since the 
rescaling factor is now 6 = 3. The transverse fiuctuations 



scale similarly to Eq. ( 4.15 ) 



C. Surface critical behavior 

At the surface we define a temperaturelike parameter 
/is = hi/ J with scaling dimension yt^ and introduce as 
above an auxiliary variable 9 to take into account the 
asymmetry of the renor mali zed couplings. 



A comparison of Eqs. {Gi) with Eqs. (C5) leads to the 
renormalized parameters 



12 2 



e / c{iis/iiy 

c — e 



(6.9) 



With the bulk Ising-fixed-point values given in 
Eqs. (6.5) and (|6^), two surface fixed points are obtained 
with 



^Ms 



= (2 + i?)-i < 1 , 
= (2 + i?) > 1 , 



,,*2 _ ,,*2 _ n 

Ms - M - R , 



,*2 



(6.10) 



At the stable fixed point = i?, with 6 = 3, the scal- 
ing dimension of the surface temperaturelike parameter 
reads 



ln(2 + R) 

hTs 



(6.11) 



a result previously-conjectured on the basis of a finite- 
size scaling study.O The same quantity at the unstable 
fixed point is given by 



ln(2 + R) 
2 ln3 



(6.12) 



The surface magnetization exponent follows as 
above from the behavior under renormalization of the 
odd components V (pk+ 1) , V ( Qk-\-2>) , and F(6fc-t-5) which 
follows from Eq s. fel| b)-(|Cl|e) and (jej). At the Ising 
fixed point Eq. (2.4) leads to 



$*(3fc+l) = $'*(/c + l) , 
$*(3fc + 2) = -i?i/2$'*(/c 1) 
$*(3fc + 3) = $'*(/c 



1) 



(6.13) 



Making use of Eq. (4.22) with 6 = 3, the scaling dimension 
of the surface magnetization is given by 



ln(2 + R) 
21n3 



(6.14) 



in agreement pith a direct calculation of the surface 
magnetization B3 Aga in a;„j, satisfies the scaling rela- 
tions ( |4.25| ) and ( |4.26| ). 



VII. PAPER-FOLDING SEQUENCE 

A. Definition and general properties 

The paper-folding sequence0 results from the recur- 
rent folding of a sheet of paper onto itself, right over 
left. After unfolding, one obtains a succession of up and 
down folds to which one associates a digit, and 1, re- 
spectively. After four steps, this process leads to the 
following sequence: 



001001100011011 



(7.1) 



The sequence on the right of the central fold is the mir- 
ror image of the left part, with each digit fk replaced 
by its complement l — fk- As a consequence, the asymp- 
totic density is poo = 1/2 and the Ising critical coupling 
is Ac = i?-i/^ 

The same sequence can be generated using the four 
letter substitutions A S{A) = AC, B S{B) = DB, 
C S{C) = DC, and D S{D) = AB with the identifi- 
cation A = 00, B=ll, C= 10, and D = 01. The leading 
eigenvalues of the substitution matrix, Qi — 2 and = 1, 
lead to a vanishing wandering exponent, a; = 0. 
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The even terms f2k, underlined in Eq. (7.1), reproduce 
the sequence itself whereas odd terms are alternatively 
and 1. Thus one has 

f2k^fk, /2fc+l = ^[l + (-l)1. (7.2) 



B. Bulk critical behavior 

The renormalization of the paper-folding problem is 
slightly more involved than the preceding ones. In the 
decimation process, as shown in Fig. ^ one eliminates 
blocks of two sites which interact alternatively via J or 
RJ. As a consequence, alternating transverse fields h'^ 
and h'p are generated at odd and even lattice sites, re- 
spectively. Furthermore, some auxiliary asymmetry pa- 
rameters are needed to retrieve eigenvalue equations with 
their original form after renormalization. Altogether the 
exact RG transformation involves six variables. 



C. Surface critical behavior 



RJ 

haJhp I 



I x \ - 



1 



FIG. 6. As in Fig. | for the paper-folding sequence. 



Using Eqs. (^) and (|D5D together with Eq. (p.4|), 
it may be verified that the normalization of here 
involves two components, $'*(2A: + 1) and ^'*{2k + 2). 
Thus the renormalization of the surface magnetization, 
based on the renormalization of the eigenvectors, be- 
comes equivalent to a direct calculation of nis- In this 
case it is more convenient to introduce, besides the sur- 
face temperature, a surface field conjugated to TOj in 
the original Hamiltonian and to study its scaling behav- 
ior. 

This can be achieved, while keeping the free-fermion 
character of the Hamiltonian, thro ugh the addition of 



a surface term —^hsCrQaf in Eq. (LS). Since there is 



no transverse field acting on the first site, ctq, which 
commutes with H, is conserved. The eigenstates of the 
Hamiltonian then belong to one of the two sectors corre- 
sponding to the eigenvalues ±1 of ctq. Thus the supple- 
mentary term takes the form =pi/iscrf and corresponds 
to a surface field ±hs acting on af, the sign depending 
on the sector. 

The decimation of the surface block described in Ap- 
pendix ^ gives the rcnormalizcd parameters 



The renormalized parameters, following from the dec- 
imation process detailed in Appendix are given by 



A2 - 1 



A2 -i?2 



/2 2 Co 

K = K — 

a 



'fi 



_ da 
Co 

di = Ki - ■ 



2 



dp 

- A2 - 1 

The fixed point values of interest for the IM are 

A* = 0, ^il = -R, A*^ = -l, 



R\ 



"/3 



R- 



(7.4) 



A linearization of the RG transformation, Eqs. (7.3), near 
this fixed point gives 



OA' 
dA 



--[{l + R){l + R-^)Y/\ 



dn0 



-R, 



-R-K (7.5) 



(7.6) 



The first line leads to the anisotropy exponent 
_ \n{R^^^ + i?-i/2) 
^" h^2 ' 

previously conjectured in Ref. whereas the leading 
eigenvalue in the linearized transformation of the tem- 
peraturelike variables, which is equal to 2, gives the cor- 
relation length exponent v = l. 



da^'jj 



ds(A2-l)2 



dsHf} 



(7.7) 



where C,s = hs/J and iis = hi/J are reduced surface vari- 
ables whereas 6 takes into account the asymmetry intro- 
duced by iig. Ci and di are the bulk parameters defined 
previously in Eqs. (|7.3| ) and ds = Ka — fi'^/ {A'^ — I)- 

Let us first consider the scaling behavior of fig, i.e., of 
the surface the rma l perturbation. With the bulk values 
given in Eqs. (7.4) one obtains two Ising surface fixed 
points with 



= {i + Ry^ < 1 



(i + i?) > 1 



,*2 



(7.8) 



The stable fixed point corresponds to fi^^ — R^ and, with 
6 = 2, the scaling dimension of the surface temperature is 
given by 



ln(l + R) 



(7.9) 



as expected from numerical results^ whereas 
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l n(l + R) 
2 In 2 



(7.10) 



at the unstable fixed point. 

The stable fixed point values of the parameters in the 
equation for the surface field variable Cs l^ad to the trans- 
formation 



(7.11) 



Thus, in the extended parameter space, there is a flow 
from C*'^ = to C*^ = +00 and the critical behavior is 
governed by the fixed point with a vanishing surfac e field 
which is unstable in the direction of . Then Eq. (7.11) 
gives the scaling dimension of the surface field for the 
Ising problem. 



ln(l + R- 
2 In 2 



or, using Eq. (7.6) 



X'rr 



ln(l + R) 
2 In 2 



(7.12) 



(7.13) 



in agreement with the scaling relations ( [4.25| ) and ( 4.26| ) 
and the analytical result of Ref. p3. 



B. Bulk critical behavior 

The quantum chain is assumed to st art o n k — 1, i.e., we 
ignore the first digit in the sequence ( ^.l| ). As indicated 
in Fig. 1^ in the renormalization proces s, o dd interaction 
terms J2fe+i which, according to Eqs. ( |8.2| ), are equal to 
J in the bulk, are eliminated so that 6 = 2. With the 
same notation as before for the reduced parameters, the 
renormalized variables follow from Eqs. (El) and (E5) 
with 



A' = ^(A^-A^-l) , A' = (A2-A2), 
A 

V\2k)^V{Ak), y'(2fc + 1) = y(4fc + 1) . 
Near the Ising fixed point, corresponding to 

A*=0, A*=-l, 
the eigenvalues of the linearized transformation 



dK' 
9A 



9A 



= 2 



(8.3a) 
(8.3b) 

(8.4) 
(8.5) 



lead to the unperturbed Ising values for the anisotropy 
and correlation length exponents, z — v — 1. 



VIII. FREDHOLM SEQUENCE 



RJ 



A. Definition and general properties 

The Fredholm sequencJ^ is obtained via substitution 
on three letters, A S{A) = AB, B S{B) = BC, and 
C ^ S{C) = CC . We start the substitution process with 
A and here, for convenience, we number the sequence 
starting on fc = 0. With A = 0, B = 1, and C = 0, after 
four iterations one obtains the sequence 



0110100010000000 



.1) 



which is the characteristic sequence of the powers of 2, 
fk being equal to 1 for k = 2P. Even underlined terms re- 
produce the sequence and odd terms, except /i, vanish. 
This gives the relations 



f2k — fk 



W-O (fc>0) 



= (8.2) 



The leading eigenvalues of the substitution matrix are 
rii = 2 and Q2 = hence the wandering exponent once 
more vanishes. 

The number of digits equal to 1, n^, grows logarithmi- 
cally with the length L, and thus the asymptotic density 
is Poo = 0. The Ising critical coupling in Eq. ( |3.9| ) keeps 
its unperturbed value Ac = 1. This aperiodic perturba- 
tion modifies the surface critical behavior but does not 
change the bulk properties, except near line defects which 
introduce local marginal perturbations in the 2D IM.E3 



X X 



xjjx 



=h I 



FIG. 7. As in Fig. ^ for the Fredholm sequence. 



C. Surface critical behavior 

With the Fredholm sequence, the decimation of the 
surface block introduces a multiplicative renormalization 
of the first component of the eigenvector such that (see 
Appendix ^) 



v'{i)^ev{i) , 



(A^-A^)(acA^-A^.R^-^^) 

(A2-A2-1)(KA2-A2i?2) 

(8.6) 



The transformation of the other parameters follow from 
Eqs. (p|) and (p|) as 



t? = tl 



(A2-A2-1)(A2-A2)i?2 



(kA2 - A2i?2)(KA2 - A2i?2 _ t2) 
(A2- A2)(KA2-A2j;2_^) 
(A2-A2-1)(KA2-A2i?2) ■ 



■7) 



As before ts = hi/h is the surface temperature and n an 
auxiliary variable generated by the transformation. 
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For the Ising fixed-point values of the bulk parame- 



ters given in Eqs. (8.4) and {kA'^)* = 0, one obtains two 
surface fixed points with 





* 2 


j.*2 
'-s 


dtl 








.*2 




2 ' 



0, 



1 

2 ' 



1 

i?2 



The first fixed point is stable when R> Rc — \/2 and, 
since t* = corresponds to a vanishing transverse field on 
the first spin, the surface is ordered at the critical point. 
The second fixed point only exists in the regime R<Rc 
where it is stable. With 6 = 2, the scaling dimension of 
the surface temperature is given by 



2-hr2' 



Ini? r- 



(8.9) 



These expressions were conjectured in Ref. ^ on the ba- 
sis of a finite-size scaling study. 

The scaling dimension of the surface magnetization 
follows from the transformation of the odd components 
of the eigenvector at the ap propriate Isi ng su rface fixed 
point. In the bulk Eqs. (^.4|), (|e|), and ( 8.3b ) lead to 



$*(2fc + 1) = $*(2fc + 2) = $'*(fc + 1) , 



fc > . 
(8.10) 



In the surface block, using Eqs. (EE) and (8.6), one ob- 
tains 



**(1) = 



$'*(!) 



$*(2) = 



t*$'*(l) 

e*R 



.11) 



Thus the normalization of gives 

£t<J'-(2fc + 0=(l + |)^ + 

oo 

+2^$'*2(fc + l) = l. 



fe=0 1 = 1 



.12) 



k=l 



According to Eqs. (B^), the coefficient of $'*2(1) is equal 
to 2 at both fixed points so that J2kLo *'*^ (^+1) = ^ - Th e 
surface magnetization transforms according to Eq. (4.23), 
i.e., like 



V2e* 



.13) 



With b — 2, this leads to the following scaling dimensions 
in the two regimes: 



a;™, = , R> V2 , 

"^'^ = 2 ~ hTI ' ^ < ^2 ' 



.14) 



as given by |-a direct calculation of the surface 
magnetizationEj The value Xm^ — when R> R^ is con- 
sistent with the vanishing surface transverse field at the 
fixed point. There is surface order when the critical point 
is approached from the low-temperature phase and, since 
the surface is one dimensional, the local magnetization 
vanishes discontinuously when the bulk disorders. When 
R<Rc the strength of the perturbed couplings is not suf- 
ficient to maintain the surface order and the transit ion i s 
continuous. In this latter case the scaling relations ( 4.25| ) 
and (4.26) are still verified. 



D. Aperiodic perturbation in the bulk 

Let us consider the aperiodic perturbation which fol- 
lows from the junction of the Fredholm perturbation in 
one half-space to its symmetric counterpart in the other, 
i.e., using the symmetrized sequence 



1000101111010001' 



5.15) 



The second half of the sequence is assumed to s tart on 
A: = l, leaving out the term A: = in the sequence (8.1) as 
for the surface perturbation. In this way one obtains a 
symmetric defect in the bulk with a vanishing asymptotic 
density so that Ac remains equal to 1. 

The simple relation between the local magnetization 
and the components of the eigenvector corresponding to 
the lowest excitation no longer holds in the bulk and one 
cannot introduce a local field term conjugated to cr^ in 
the Hamiltonian (1.2) without breaking its free-fermionic 



character. Thus we shall only consider the renormaliza- 
tion of the local temperaturelike variable Td = hi/h. For 
the other p arameters we keep the same notation as in 
Sec. VIII B for the surface. 



RJ 
h 1 J 

xl 1 x|x 
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X X 


X 


x|x' 
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FIG. 8. As in Fig 


□ for the Fredholm defect in 


the bulk. 



The decimation of the central block of eigenvalue equa- 
tions is illustrated in Fi g. ^ . The rcn ormalized local vari- 
ables follow from Eqs. (E7) and (ES) as 



Ai?2(A2~A2) 



(kA2 - A2i?2)2 
A2 - A2 

1 



A2-A2 



1 



K(^^A2-A2i?2^ 



(kA2 - A2i?2) 



-JA2 



(8.16) 



In the critical surface, with A* = and A* = — 1 at the 
Ising fixed point, the RG transformation of the local vari- 
ables takes the form 
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Td 
i?2 



1 



i?4 



.17) 



When i? = 1, T^j^ = Td, which leads to a hne of fixed 
points parametrized by Td, k* — 1 + Td. The scaUng di- 
mension of the local temperature vanishes as expected 
for a thermal line defect in the 2D IM.L3 

When i? 7^ 1 , two fixed points are obtained with 



dTd 



1 



= 0, 



i?2 



2i?2 - 1 



0. (8.18) 



In the critical surface, the fixed point at = is sta- 
ble when R>1. The transverse field at the center of the 
defect vanishes at this fixed point. Thus the defect is 
ordered at the critical point, like for the surface, but the 
critical value of R is now i?c = 1 instead of ^/2. There 
is no need to compensate for missing bonds as it is the 
case at a surface. When R< 1, the second fixed point 
at — +O0 becomes stable and leads to a second-order 
transition at the defect. The appropriate scaling field td, 
associated with the defect temperature, which vanishes at 
the fixed point is now tJ^. Hence, with b 



.1? 



= 2, Eqs 

give the following scaling dimensions in the regimes of 
first- and second-order local transition, respectively 



InR 

yt, = -2 , R>1 

^ " In 2 ' 

Ini? 

Vt, = 2 -r-r , R<1 . 
ytd In 2 



(8.19) 



IX. 



RELATIONS BETWEEN ISING MODEL 
CRITICAL EXPONENTS 



Apart from the correlation length exponent = 1, all 
the critical exponents obtained for the different aperiodic 
models are varying with the amplitude of the modulation, 
and thus the critical behavior of these models is nonuni- 
versal. However, some kind of "weak universality" still 
holds and there are relations between critical exponents 
which follow from the fact that the_systems at the criti- 
cal point obey anisotropic scaling.o A detailed analysis 
of the scaling behavior can be found in Ref. |3^. 

One can notice other exponent relations which are spe- 
cific for the marginally aperiodic IM's. One such relation 



connects the scaling dimensions of the energy and mag- 
netization densities at the surface as 



Z I 2i3j 



(9.1) 



It follows from Eq. ( 4.25 ) and anisotropic scaling. It 
was conjectured in Ref. |33| on the basis of an assump- 
tion for the scaling behavior of $(1) for low- lying excita- 
tions. One can find another relation which surprisingly 
connects bulk and surface quantities in the form 

where Xm, is the scaling dimension of the surface mag- 
netization on the right-hand side (RHS) of the system. 
Here we argue that the relation in Eq. (9^) is generally 
true for marginally aperiodic layered IM's. 

In the following derivation, we consider the quantum 
Ising chain Hamiltonian TL given in Eq. (1.2) with homo- 
geneous transverse field /ifc = 1. The dynamical exponent 
of the model z is related to the scaling behavior of the 
lowest gap of the spectrum of the critical Hamiltonian in 
the form 



Ei-Eo^ Ai 



(9.3) 



in a finite system of size L. 

The asymptotic size dependence of Ai [L) is calculated 
in the following approximation. First we determine the 
leading k d epe ndence of the eigenvectors $i (fc) and ^'i (A:) 
from Eqs. (2.3) in such a way that the RHS's of the equa- 
tions are omitted. This approximation is justified, at the 
critical point or in the ordered phase, by the fact that the 
second difference operators on the LHS of the equations 
are 0{L~^) whereas on the RHS is 0{L^^^) with z> 1 
for marginal aperiodic systems at criticality or exponen- 
tially small in the ordered phase. In this approximation 
we obtain 



fe-i 



$i(L + l-fc)~$i(L) J](-Ai_, 



2 



i=\ j=l 



fe-1 



*i(A:)~*i(l)[](-A, 



1=1 



(9.4) 



i=i j=i 

Then the size dependence of Ai is estimated from the 



linear equations in Eqs. (2.2) as: 



ML) 
Ai(L) 



^i(l) 
$1(1) 

^i(^) 
«'i(L) 



L-l 



1=1 

L-l 



L-l i 



^ ' 2 — 1 



1=1 j=i 

L-l i 

i=i j=i 



-2 



(9.5a) 
(9.5b) 
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Multiplying both sides of Eqs. (9.5a) and (9.5b), one ar- 
rives at the result 



L-1 



Ai(L) ~ m,(L)m,(L) []A-i 



(9.6) 



1=1 



where the finite-size surface-tnagnetizations on both sides 
of the system are given byila 



ms{L) 
rns{L) 



L-1 i 

i=i j=i 

L-1 i 

i=i j=i 



-1/2 



-1/2 



(9.7) 



The relation in Eq. (9^), which connects the asymptotic 
behavior of the lowest excitation energy and the finite- 
size behavior of the surface magnetizations, is valid for 
general distribution of the couplings, provided the lowest 
gap in the system goes to zero faster tha n 1/L. 

In the following, we apply Eq. ( |9.6| ) to marginally 
aperiodic systems at t 
ing to rigorous results. 



critical point where, accord- 
nfji^(^»)c = 0(1) for aperi- 
odic perturbations leading to a shift of the critical cou- 
pling. The finite-size surface magnetizations behave as 
rUsciL) ~ L~^™» and rnsc{L) ~ and thus, from 



Eqs. ( |9.6D a nd ( |9.3D , one obtains the scaling relation given 
in Eq!7^). 

The aperiodic sequences studied in this paper which 
change the bulk critical behavior are of two kinds: ei- 
ther symmetric with Afe = Xl-u (period doubling) or 
such that a perturbed coupling at k corresponds to an 
unperturbed coupling aXL — k, which leads to p^o = 1/2 
and, according to Eq. ( ^j| ), [Xk{R)]c^[\L-k{R~^)]c (pa- 
per folding, three folding). For symmetric sequences, 
rnsc{L, R) = msc{L, R), and therefore Xm, = Xm,- Oth- 
erwise, rfisc{L,R) = msc{L, R~^) and, consequently, 
XnisiR) — XmS^~^)- Thus the knowledge of a single 
exponent Xm^ (R) is sufficient to obtain all the varying 
exponents studied in this paper. Furthermore, for the 
period-doubling sequence x„i^ (R) is symmetric under the 
exchange of R into R~^ according to Eq. (4.24). It fol- 
lows that for all the aperiodic sequences one may rewrite 
the nonuniversal anisotropy exponent in Eq. (9.2) as 

z{R)=X,nAR)+^mAR~')- 

For marginal aperiodic sequences which do not change 
the bulk critical behavior, i.e., leave z — I and Ar = 1, 



the scaling relation (9.1) does not hold. In this case 
]^^^^(Ai)c = -R~"^ with the number of perturbed cou- 
plings growing logarithmically with L. For the Fred- 
holm sequence = In L/ In 2 so that the product of 
the couphngs in Eq. ( |9.6[ ) scales as ^v"'"-'^/'"^. When 
the left surface is ordered at the critical point, i.e., for 
R > Rc = v^, and the right surface is free, we have 
Xm^ — 0, Xm^ —1/2 and the lowest excitation does not 



scale as L -'^ like the rest of the spectrum, but with a 
continuously varying exponent: 



Ai ^ 2^-l/2-ln-R/ln2 



(9.8) 



When the surface magnetization vanishes at the c ritica l 
point {R<Rc), the R dependence of Xm^ in Eqs. ( 8.14 ) 
just compensates that appearing in the product of the 
couplings and one recovers the normal behavior for 
Ai. 

Finally, according to Eq. ( 4.26| ), the scaling dimension 
of the surface magnetization Xm^ is equal to the scaling 
dimension of the surface temperature y^^ at the unstable 
fixed point. This last relation is a cnaaequence of the 
self-duality of the Ising quantum chainJl3 Using the dual 
Pauli spin matrices defined through 



~k~k+l > 



(9.9) 



the original Hamiltonian in Eq. (1.2) is transformed into 
its dual: 



L-1 



n 



'fc-l-1 ) 



(9.10) 



fc=i 



k=l 



with a vanishing transverse field on the fir st spin. As al- 



ready shown at the begining of Sec. VII C, rf commutes 



with 7i and may be replaced by its eigenvalues ±1. Thus 
in the surface term — i/iirf rf = ^F^/iirf , hi = htg now 
plays the role of a surface field acting on rf . 

The unstable fixed point at t* = 0, with its associated 
scaling dimension yt^ , governs the critical behavior of the 
dual surface magnetization. In the duality transforma- 
tion the couplings Xk^h/ Jk are changed into A^^ so that 



the surface magnetizations on both sides of Eqs. (9.7) are 



exchanged. It follows that the scaling dimension of (r| 
is Xm^ and the dimension of the surface field is given 
by — z — Xm^ , which, according to Eq. ( |9.2|) , leads to 
Eq. (^^261). 



X. CONCLUSION 

In this paper we have presented a unified statistical- 
mechanical description of the IM and the DW on layered 
two-dimensional lattices, taking the extreme anisotropic 
limit for the IM. The critical properties of the two prob- 
lems were deduced from the scaling behavior of the spec- 
trum of the transfer matrix of the DW, which is studied 
through exact RG transformations. For a given value 
of the aperiodicity parameter R the RG transformations 
have two nontrivial fixed points, as shown in Fig. ^. The 
bottom of the spectrum scales to the IM fixed point, 
which controls the critical behavior of the IM, whereas 
the top of the spectrum scales to another fixed point, 
which describes the critical properties of the DW. 

The aperiodic sequences we considered have different 
effects on the critical properties of the two models accord- 
ing to Luck's relevance-irrelevance criterion described in 
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Sec. [II 



^ For the IM, the crossover exponent in Eq. (3.7) 
1 is 0iM = O whereas for the DW, with z/_l = 1/2, 



with V 

it is (/)Dw = 1/2. Consequently, the nonperiodic pertur- 
bation is marginal at the homogeneous {R~l) IM fixed 
point whereas it is relevant for the i? = 1 DW fixed point. 
These statements are in accordance with the exact re- 
sults. 

For the IM, the marginal perturbation creates a line of 
fixed points, which is parametrized by i?, and the crit- 
ical properties are continuously varying, even at i? = 1 . 
The nonperiodicity also induces a continuously varying 
anisotropic scaling behavior. However, the different vary- 
ing exponents are not independent: Knowledge of the 
scaling dimensions of the surface magnetization is suf- 
ficient to completely describe the nonuniversal critical 
behavior studied in this work. 

Considering the DW problem, here the line of fixed 
points is discontinuous at i? = 1 , in accordance with the 
relevant nature of the perturbation. For the hierarchical 
models the line of fixed points is characterized by finite 
coordinates and the corresponding critical behavior is of 
power-law form with i?-dependent exponents. On the 
other hand for the aperiodic models (period doubling, 
three folding and paper folding) the line of DW fixed 
points is shifted to infinity and the scaling behavior is 
anomalous: The transverse fluctuations of the walk grow 
on a logarithmic scale. 

Finally we discuss the local critical behavior at ex- 
tended defects, located either at the surface or in the 
bulk, which are generated by the Fredholm sequence. In 
both cases two fixed points exchange there stability at a 
critical value Rc of the modulation amplitude R. This 
critical value separates two regimes for the local transi- 
tion: For R> Rc the local magnetization vanishes discon- 
tinuously at the bulk critical point while for R< Rc the 
transition is continuous. In both cases one obtains criti- 
cal exponents which vary continuously with the marginal 
parameter R. _ 

It has been already noticecCj that the surface Fred- 
holm perturbation is_closely connected to the Hilhorst- 
van Leeuwen modelO In the same way, the Jaiilk Fred- 
holm defect is connected to the Bariev modeLEa In these 
models, the perturbation of the couplings decays as a 
power of the distance I from the center of the defect with 
6X{l)=al~^ in the marginal case, for the 2D IM. 

The varying exponents obtained analytically and nu- 
merically in Ref. ^ for the surface Fredholm perturba- 
tion as well as th ose otained via exact RG transforma- 
tions in Sec. VIIl for the surface and bulk Fredholm de- 



fects can be put in correspondance with the exponents of 
the Hilhorst-van Leeuwen and Bariev models with a re- 
placed by In i?/ In 2. Up to now, the values of the Bariev 
model's exponents had been conjectured on the basis of 
conformal methods using gap-exponent relations after a 
conformal transformation of the inhomogeneous infinite 
system onto an inhomogeaeous infinite strip with peri- 
odic boundary conditions.^ Our RG results for the bulk 
Fredholm defect and the correspondence between both 



models strongly support the validity of this procedure. 
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APPENDIX A: RENORMALIZATION OF THE 
PERIOD-DOUBLING SEQUENCE 

The renormalization of the period-doubling sequence 
in the bulk involves the following set of equations: 

XR^^" V{8k) - A U(8A:+1) + V{8k+2) = , (Ala) 

V{8k+ 1) - A U(8fc+2) + A U(8fc+3) = , (Alb) 

A U(8fc-|-2) - A U(8fc + 3) + U(8fc+4) = , (Ale) 

U(8fc + 3) - A U(8fc+4) + XR Vi8k+5) = , (Aid) 

Ai?U(8fc + 4) - AU(8fc+5) + U(8fc+6) = , (Ale) 

U(8fc + 5) - A U(8fc-|-6) + A U(8fc+7) = , (Alf) 

A U(8fc-|-6) - A U(8fc + 7) + A U(8fc + 8) = , (Alg) 
U(8fc + 7) - A U(8fc-^8) + XRf^x^ U(8fc+9) = . (Alh) 

Among these eight equations we eliminate the six central 
ones, which amounts to rescale the system by a factor 
b — A. This is accomplished by evaluating U(8fc-|-2) and 
V{8k + 7) as functions of U(8fc-|-1) and V{8k + 8), in the 
linear system given by Eqs. (|Al|b)"( Alg), with the result 



V{8k+2) = —V{8k+l) 
c 



RX^ 



V{8k+7) 

c 

c = A2(A2-A2 
d^{A^~ 1)(A2 



U(8fc+1) 



RX^ 

c 

dl 



c 



V{8k + 8) , 
U(8fc + 8) , 



2\2 



If ~X'R'{K' ~X') 
X^ -l)-X^R^(k^ -X^) 



(A2) 



Inserting these values into Eqs. (Ala) and ( Alh| ), after 
multiplication by c/{RX'^) we obtain 



^ Rf^" V{8k) - A ^ 



i?A2 
V{8k 



i?A3 



U(8fc + l) + U(8fc-l-8) = 0, 



-D- A^^TT^ U(8fc+8) 



i?A3 



+ V{8k+9) = , (A3) 



which are the renormalized equations 
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X'Rf" V'{2k) - A' V'{2k + 1) + V'{2k+2) = , 
y'(2fc+l)-AV(2fc+2) + A'i?-^'^+iy'(2fc+3) =0, (A4) 



after rescaling by 6 = 4. Noticing that, according to 
Eqs. ( |4.2D , fik = 1 — /2fe = fk, R remains unchanged and 
one obtains the RG transformation as given in Eq. (4.3). 

At the surface, in terms of the redu ced variables, the 
same set of equations as in Eqs. (Al) with fc = is ob- 
tained, except for the two first equations which now read 



ts 



-AV{1) +0tsV{2) = 0, 
V{1) -AV{2) + XV{3) = 



(A5a) 
(A5b) 



The auxihary variable 9 is needed to take into account 
the asymmetry resulting from the renormalization after 
one step. In this way the variables A, A, tg, and build 
a closed set under renormalization. 

As above, the components V{2) and V{7) can be de- 
duced from the six central equations and read 



V{2) 
Vi7) 



dAts 
RXH 



-v{i) 



RX^ 
c 

dA 
f — 

c 



V{8), 
ViS), 



(A6) 



wher e c an d d are defined in Eqs. ([A2|). Equations (A5a) 
and ( Alh| ), after multiplication by appropriate factors, 
then give 



dt2 



■vil) 



F(8) = , 



V{9) = 



(A7) 



These equa tions give the renormalized forms of Eqs. 
( A5a ) a nd (A5b ) and provide the RG recursions given 
in Eqs. (O). 



APPENDIX B: RENORMALIZATION OF THE 
HIERARCHICAL SEQUENCE 

In the bulk, the set of eigenvalue equations we consider 
is the following: 

Ai?" V^(2m") - AF(2m" 4- 1) -I- V{2m'' + 2) = , (Bla) 
y(2m"-|-l)-AF(2m"-|-2) + AF(2m"-|-3) = 0, (Bib) 
AP"(2m" -I- 2) - Ay (2m" + 3) + V(2m" +4) = , (Blc) 

AF(2TO" + 2TO-2)-Ay(2m"-h2TO-l) + 

+y(2m" + 2m) = 0, (Bid) 
y(2m"-F2m-l)-A\/(2m" + 2m) + 

+XR \/(2m" + 2m + 1) = . (Ble) 



Among the 2m equations one eliminates the 2to — 2 cen- 
tral ones, which amounts to rescale the system by a factor 
of b — m. Then we are left with two equations between 
the components y(2m"), V{2m'^ + 1), V{2m'^ + 2m) and 
y(2m" + 2m-Hl) of the form 



Ai?" 



V{2m" 



A-s 



y(2TO" + l) + F(2m" + 2m) = 0, 



T/(2m" + l) - ^^y(2m" + 2TO) + 
r 

+ hllv(2m'' + 2m + 1) = . (B2) 
r 

Here r = A""V-D2m-2, whereas s = -D2m-3/ D2m-2, 
and I?2m-2 denotes the (2m — 2) x (2m— 2) determinant 



D 



2m-2 



-A A 
A -A 1 

1 -A A 



1 -A A 

A -A 



(B3) 



while D2m-3 is the lower central minor of D2m-2- Then 
from Eq s. (B^ ) we deduce the RG transformation given 
in Eqs. (js!?). 



APPENDIX C: RENORMALIZATION OF THE 
THREE-FOLDING SEQUENCE 



The eigenvalue equations take the following form in the 
bulk: 

Rf-^>'V{<6k)-KAV{<6k+l) + ^V{Qk + 2) = Q, (Cla) 

^iV{&k + l)--V{&k + 2) + V{Qk+?,) = {) , (Clb) 

V{Qk + 2)~KAV{Qk+i) + ^iV{Qk+A):^Q, (Clc) 

yy{&k+i)--V{&k+A) + RV{<ok+^) = {), (Cld) 

RV{Qk+A)-KAV{<ok+b) + ^lV{<ok+Q) = Q , (Cle) 

A f 

/iy(6fc+5) t/(6fc+6)-fi?^^'=+^ y(6fc+7)==0. (Clf) 



Equations (Clb)~(Cl3) can be used to write 



X^(6fc+2) 



ckA 
e/i 



y(6fc-f 1) 



V{Qk+^) = ^V{Qk + l) 



R^i^ 

e 

dA 



e 

,2 d2\ 



.2/T2 



V{6k + 6). 
V{6k + 6) , 



c = ^i^A^-^l^-R^), d = fi\A^-fx^-l). 
e=(A^-l)(A^-R^)-ii^A\ 



(C2) 
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which, inserted into Eqs. (Cla) and (Clf), lead to the 
renormaUzed equations 

Rf^^^V(6k)-KA(l--) y(6fc+l) + ^V(6fc+6) = 0, 
V e/ e 

^V{6k + 1)-- ( 1--] V{6k+6) + 
e K \ e J 

+Rf'"'+Wi6k + 7) =0. (C3) 



Since f^k = fk according to Eqs. (6.2), these equations 
take the form 

Rf'-V' {2k)- k'A'V' {2k+l) + fi'V' {2k+2) ^0 , 

n'V'{2k + l)-^V'{2k+2) + R^'''+'V'{2k+3) ^0, (C4) 



with the renormahzed vari ables g iven in Eqs. ( p. 3D 

At the surface, Eqs. (|Cllc)-(|Cl|') with fc = have to be 
supplemented by 



-KAV{l) + dnsV{2) = , 
^^(l)_^y(2) + y(3) = 0. 



(C5) 



Rewriting V{2) and V{5) as functions of 1^(1) and V{6) 
in the first and last equations of the surface block, one 
obtains the renormahzed equations 



^(6) = 0, 



^ni)-^(i-f)n6)+n7)=o. 



(C6) 



A comparison with Eqs. (C5) leads to the renormahzed 
parameters given in Eqs. (6.9). 



when the central interaction is RJ. /Zq = ha/ J and 
/i_a — hp/ J are reduced temperaturelike parameters and 
A the reduced eigenvalue defined before. Since in each 
block two sites out of four are eliminated, lengths are 
rescaled by a factor 6 = 2. 



The two intermediate equations in Eqs. ( pi| ) and ( |D2| ) 



give 



V{%k+2) = ^^y(8fc + l) + J^y(8fc + 4) , 
y(8A;+3) = J^V^(8fc + l) + ^C^y(8A;+4) , 

A2 - 1 A2 - 1 

y(8A:+6) = .^^l^(8fc+5) + J^y(8fc+8) 



A2-i?2 



A2-i?2 



y(8fc+7) = J^^n8fc+5)+ ~ '^^'^ y(8fc+8), (D3) 



A2-i?2 



A2-i?2 



whic h ca n be used in the first and last lines of Eqs. ([D1| ) 
and (D2), together with the first relation in Eqs. ([7.2|), 
to write the renormahzed equations 



Rf V' (4fc) - k'^A'V' (4fc+ 1 ) + /x^T/' (4fc+ 2) = 0, 
A' 

K 



li'aV {Ak+1) -F'(4fc+2) + i?/^'=+iV"'(4/c+3) = 0, 

i?^=''-+iy'(4/c+2)-K'A'y'(4fc+3)+/i^^^'(4fc+4) = 0, 



n'pA' 



ti'pV\^k+i) ^T/'(4fc + 4) + i?/^'=+2y(4A:+5) = 0. (D4) 



Here the components of the eigenvectors transform ac- 
cording to 



APPENDIX D: RENORMALIZATION OF THE 
PAPER-FOLDING SEQUENCE 

The following blocks have to be considered, 
i?/«T/(8fc)-K„KAF(8fc + l)+AiaV^(8fc+2) = 0, (Dla) 

/i„F(8fc+l)--y(8/c+2) + V(8/c+3) = 0, (Dlb) 

K 

y(8fc+2)-KAy(8A;+3) + /i/3V^(8fc+4) = 0, (Die) 

/i/3y(8A:+3)-^y(8/s+4) + i?^^'=+^F(8fc + 5) = 0, (Did) 

when the central interaction is J and 

i?-^"+2F(8fc + 4)-K„KAF(8fc + 5)+/i„V^(8fc+6) = 0, (D2a) 



AiaV^(8fc+5)--y(8/c+6) + W(8fc + 7) = , 



V'{Ak) = V{Sk) , V'{Ak+l) = V{Sk + l) , 
y'(4/c+2) = y(8fc + 4), y'(4/c+3) = y(8fc + 5). (D5) 



The renormahzed parameters are given in Eqs. ( |7.3| ). 

The surface field h^ = JCs introduces a supplementary 
equation in the surface block which now begins with 



Kf3A 



V{0) + CsVil) = , 



Cs v{o) - KaKAv{i) + efj.sV{2) = 

^ni)--n2) + n3) = 0, 

K 



(D6) 



(D2b) 



RVi8k+6)-KAV{8k+7) + ^i0Vi8k+8)=O, (D2c) 
fif3V{Sk+7)-'^V{8k+8) + Rf-"'+''V{8k + 9) = 0, (D2d) 



where, as before, fig = hi/ J is a temperaturelike surface 
variable and an auxiliary variable. 

The first equation in Eqs. ( p^ ) gives the value of 1^(0) 
which can be used in the second equation to obtain a 
surface block in its standard form: 
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v(i) + efisV{2)^o, 



f^V{l)--V{2) + V{3)^0 

K 



(D7) 



The two remaining equations are given by (Die) 
and (Did) with k = 0. As usual, writing V{2) and 1^(3) as 
functions of ^(1) and ^(4), one obtains the renormaUzed 
equations 

'd.nl ( 1 - ^ Vil) + V{A) ^ , 



4(A2-1) 



V^X/(4) + y(5)=0 

0(A2-1) K 



(D8) 



wher e ds — Ka — A*s/(A^ ~ !)■ When compared to 
Eqs. (D7) they lead to the renormalized parameters given 

in EqsT^Ll)- 



APPENDIX E: RENORMALIZATION OF THE 
FREDHOLM SEQUENCE 

We have to consider the following block of equations: 

Ai?^^'=t/(4fc)-AP"(4fc+l) + F(4fc+2) = 0, (Ela) 

T/(4fc+l)-AV^(4fc+2) + AV^(4fc+3) = 0, (Elb) 

AF(4fc + 2)-AT/(4fc+3) + F(4fc+4) = 0, (Elc) 

F(4fc+3)-AT/(4fc+4) + Ai?-''^'=+^T/(4fc+5) = 0. (Eld) 



Equations (pi|b) and (pill|c) give the eigenvector compo- 
nents 

F(4fc+2) = ^^^^^ — V{4:k+1) + ^ — F(4fc + 4), 



A2-A2 



A2-A2 



V{'ik+3) = ^ — V^(4fc + 1)+ ^ ^ F(4fc + 4), (E2) 
A2-A2 A2-A2 

which can be used to rewrite the first and last equations 



{P - \^)Rf^^ V(Ak) - ^ (A2 - A^ - l)F(4fc+ 1) + 



+F(4fc + 4) = 0, 

V(4fc+l)-^(A2-A2-l)l^(4fc+4) + 

+ (A2-A2)i?^^'=+^F(4fc+5) = 0. (E3) 

Since /2/c = fk and f2k+2 = fk+i , the renormalized equa- 
tions keep their original for m with the transformed pa- 
rameters given by Eqs. (B.3). 
The surface block reads 



-AV{l) + tsVi2)^0 , 
t, V{1) - kA V{2) + XR V{3) ^ 
XRV{2) -AV{S) + ViA) = , 
V{3) -AV{4:) + XRV{5) = . 



(E4a) 
(E4b) 
(E4c) 
(E4d) 



where the auxiliary variable k takes into account the 
change of the intermediate inter acti on which is no w X R 
instead of A for the bulk Eqs. (^). From Eqs. (||b) 
and (E4c) we deduce the eigenvector components 



Vi2) 
Vi3) 



t.A 



kA2 - A2i?2 

tsXR 

kA2 - A2i?2 



V{1)- 
Vil) 



XR 



kA2 - A2i?2 

kA 

kA2 - A2i?2 



y(4) , (E5) 



which are used to rewrite Eqs. (E4a) and (E4d) as: 



-A' 9V{1) +t'^V{4:) = , 

t', 0V{1) - k'A' V{A) + X'RV{5) = 



(E6) 



In these renormalized equations 6 can be interpreted as 
a renormalizatio n fa ctor for V{1) which transforms ac- 
cording to Eqs. (|8.6|) . The R G e quations for the other 
parameters are given in Eqs. ( ^.7[ ). 

In the bulk, the block of eigenvalue equations to be 
renormalized, corresponding to the center of the defect, 
is the following: 



XRV{-2) - AF(-l) + V^(0) = , 
V{-1) -AV{0) + XRV{1) ^ , 
XR V{0) - kA V{1) + Td V{2) = 
Td V{1) ~ kA V{2) + XR V{3) = 
XRV{2) -AV{3) + V{4) = , 
F(3) - Al/(4) + Ai?y(5) = . 



(E7a) 
(E7b) 
(E7c) 
(E7d) 
(E7e) 
(E7f) 



With A = in Eq. (E7c), the three last equations dif- 
fer from the surface Eqs. ( |E4|) only through the auxiliary 
factor K in the first one which is necessary to preserve 
the symmetry of the block. 

Expressing T^(0) and V{3) in terms of y(— 1) and V"(4), 
the first and last equations of the block take the same 
form as the two central ones, 

X'R V{-2) - k'A' V{-1) + t'^ V{4) = , 
TdV{-l)~K'A'V{4) + X'RV{5)^0, (E8) 

with the renormalized local parameters given by 
Eqs. (|8T6|). 
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